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The stress versus strain curves in dense colloidal dispersions under start-up shear flow are investigated com- 
bining experiments on model core-shell microgels, computer simulations of hard disk mixtures, and mode 
coupling theory. In dense fluid and glassy states, the transient stresses exhibit first a linear increase with 
\ the accumulated strain, then a maximum ('stress overshoot') for strain values around 5%, before finally ap- 

proaching the stationary value, which makes up the flow curve. These phenomena arise in well-equilibrated 
systems and for homogeneous flows, indicating that they are generic phenomena of the shear-driven transient 
structural relaxation. Microscopic mode coupling theory (generalized to flowing states by integration through 
the transients) derives them from the transient stress correlations, which first exhibit a plateau (corresponding 
to the solid-like elastic shear modulus) at intermediate times, and then negative stress correlations during 
the final decay. We introduce and validate a schematic model within mode coupling theory which captures 
all of these phenomena and handily can be used to jointly analyse linear and large-amplitude moduli, flow 
curves, and stress-strain curves. This is done by introducing a new strain- and time-dependent vertex into 
the relation between the the generalized shear modulus and the transient density correlator. 

c3' 



On 



-a 

o 
o 



PACS numbers: 82.70.Dd, 83.60.Df, 64.70.Pf 
Keywords: Colloids, Nonlinear rheology, Glass transition 



I. INTRODUCTION 



" Understanding the microscopic mechanisms dominating the mechanical response of viscoelastic liq uids and disor - 
dered solids remains an open and challenging problem of fundamental nature in the material sciences (|Larsonl[l999l) . 



■ Approaching vitrification, the viscous relaxation time increases from microscopic values characteristic for dilute or 
high-temperature systems to arbitrarily large values which easily exceed the experimental measurement time. Thus, 
already very small externally-applied flow rates can compete with the intrinsic structural relaxation and can be used to 
gain deep insights. Colloids recently have moved into the focus as they make possible microscopic in vesti gations of the 
shear- driven particle rearrangements close to the glass tra n sition ; for example iBesseling et al\ ( 20071 ) and lZausch et al\ 



( 20081 ) employed confocal microscopy and iPetekidis efal\ (|2002t ) li ght scattering. Moreover, the nonlinear rh eology of 



colloidal dispersions is of interest in itself and for many technological applications ( Mewis and Wagner[|2012l) . Slightly 



polydisperse colloids, which mimic hard spheres, have repeatedly been employed as a most simple mode l system ex- 
hibitin g gl assy arrest where th e nonli near rheology can be studied in detail; see e.g. the recent works bv lPham et all 
(|2008l ) and ISiebenburger et al\ ( 20121 ) which provide access to the broader literature. 



The transient response of a viscoelastic material to shear deformation which is switched-on at time zero and then 
increases with a constant shear rate ('start-up flow') provides d etailed insight into the tim e -dependent vi scoelasticity. 
While t his is a familiar rheo logical technique in soft materials (jMewis and Wagnerl . l2012r . iLarsonl I1999I ). it was also 
used bv lHarmon et al\ ( 20071 ) to investigate metallic glasses, with broadly similar findings, raising the exciting question 



about the fundamental relations between 'hard' and 'soft' glasses. In simulations of glasses at very low temperatures, 
individual rearrangements wh ich c ause stress relaxation upon shear deformation could be identified and studied by 
Malone v and Lemaitr 1 (|2006t) and lTsamados et all (|2009f ). Intriguingly, already for small deformations (strains), the 



stress in small systems does not increase linearly but exhibits characteristic saw-teeth indicative of the local shear 
rearrangements. Their microscopic nature in disordered materials is not as well understood as in crystalline materials. 
Detailed discussions have up to now only been possible within macroscopic nonequilibrium thermodynamics; in the 
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shear transformation zone theory bv lFalk and Langerl ( 201 lh and in older flow-defect theories bv lSpaepenl ( 19771 ). At 
higher temperatures or at not too high densities in colloidal dispersions, there appears to be an initial linear regime 
i n the stress versus strain curve, whose slope measures the elastic shear modulus of the amorphous solid material 
( Pham et all [2008) . At long times after switching on the shear rate, when the strain accumulates to the order of 
100%, the s ystem approaches a s teady state where an unique, history independent stress can be observed, as verified 
in detail by ISiebenburger et ali ( 20091 ). The resulting flow curves (viz. stationary stress as function of shear rate) 
exhibit a dynamic yield stress whenever the quiescent state is glass like ( Siebenbiirger et ali . I2009L |20 12h . Inbetween 
the linear small strain and the asymptotic steady regime, the stress versus strain curve often exhibits a maximum 
(stress overshoot) which sometimes is taken as measure for a kind of 'yield stress'. Yet, while the linear response and 
stationary state of glassy materials often can be considered simple in the sense that it is homogeneous (viz. exhibits a 
spatially constant shear rate) and history independent (viz. independent of the 'waiting time' inbetween preparation 
of the state and swit ch-on of the shear), the intermediate stress-overshoot regime is more complicated. In gel-like 
colloidal dispersions, iDivoux et al. (2011) observed inhomogeneous flow (wall slip and shear banding) beyond the 
linear regime causing the stress overshoot when the materi al failed locally. Inhomogen eous flow is quite prominent in 
shear-driven aggregated dispersions as reviewed recently bv lMansard and Colinl (l2012f). and has also b een obtained in 
theore tical approaches to stress-strain curves like in the flo w-defect approac hes bv lSteif et ali (jl982al) -. [M 

anning et ali 

(2007 ) and in the mesosc opic soft glassy rheology model bv lMoorcroft et aZl(l201ll). The soft gl assy rheology approach 
bvlFielding et al. I (120001) and the nonlinear Langevin approach by IChen and Schweizerl (|2010f ) (recently reviewed by 
IChen et ali (|2010] )) also predict the second mentioned subtle property of the stress overshoot, namely its dependence 
on the time since pre paration of the gla ssy sample (viz. aging). This has been investigated in detail in a paste 
colloidal suspension bv lDerec et all dl003), and in computer simulations of a glass forming Lennard- Jones mixture by 
IVarnik et all ()2004l) . who also observed shear banding in flows with narrow rheometer gaps. 

In the present contribution, we study the stress versus strain curve of a model colloidal dispersion, of a simulated 
binary mixture of hard disks, and in microscopic mode coupling theory in order to explore the possibility for a stress 
overshoot arising generically in viscoelastic materials in homogeneous flow and independent of the age of the sam- 
ple. While stress overshoots under such simple conditions are known in polymeric materials and from a number of 
the mentioned theoretical approaches, we focus on the question whether an approach (viz. mode coupling theory) 
starting from the particle interactions and describing broadly visco-plastic phenomena, can capture it. We present 
age- independent data from experiments and simulations below and above the glass transition packing fraction of mode 
coupling theory, and compare with theory. The experimental system is a well studied slightly polydisperse core-shell 
microgel dispersion, which allows for a good control of the equilibrat ion excluding aging - effects . Its flow curves and 
linear shear moduli close to the glass transiti o n wer e determined by ISiebenbiirger et 'ali (|2009h . and its large strain 
amplitude oscillatory stresses bv lBrader et ali ( 2010l). While the homog eneity of the flow in the system could not yet 
be investigated, the confocal microscopy study by Zausch et al\ ( 2008J) established a homogeneous flow profile in a 
fluid colloidal dispersion and also showed that the stress overshoot is related to the single particles' mean squared dis- 
placements. Brownian dynamics simulations of a binary mixture of hard disks in two dime nsions are included a s there 
the flow profile can be recorded and homogeneity can be tested. Previous investigations bv lWevsser and Ha i nail ( 201 lh 
determined the glass transition density of the system, which we use to exclu de age dependences. Moreover a number o f 



rheological studies were already performed on this computer glass former ( Henrich et al. , 20091 ; iKruger et aZ.L [20TTh . 



The theoretical investigations use the first principles mode coupling theory bv lGotzd (12009) and coworker s generalized 
to flowing systems in the integration through transients approach suggested by iFuchs and Catesl (120021). Its re sults 
for two-dimensional hard disks were compared to the Brownian dynamics simulations by Henrich et all ( 20091) and 
iKriiger et al. and in the present contribution are the basis for the development of a simplified schematic 

model. Its advantage is the possibility to jointly analyze data of the linear shear moduli, flow curves and stress vs 
strain curves in a variety of glass forming fluids, whose microscopic description is yet out of reach. It generalizes a 
schematic model from the literature, which has already successfully been used to capture flow curves and linear and 
nonli n ear moduli from different glass formers, includ i ng of the experimental sy stem we investigate ( Crassous et ali 
12008k ISiebenbiirger et ali [20091 iBrader et aUl2010l : ISiebenbiirger et ali [20121 ). 

The main outcome of our investigations will be that the transient shear stress after switching on a fixed shear rate 
at time zero generically exhibits a maximum before settling on the stationary flow-curve value. It arises from negative 
stress correlatio ns that build u p wh en th e intrinsic structural relaxation gets strongly affected by the shear driving. 
As discussed bv lKriiger et ali ( 20101 ) and lZausch et ali ([2008) , the transient dynamics initially follows the quiescent 
one, while it coincides with the stationary dynamics under flow at long times. Inbetween, it exhibits super-diffusive 
m otion and negative str e ss cor relati ons, which are h e ralded by anisotropic structural correlations seen in simulations 
by IZausch and Horbachl ( 20091 ) and iKoumakis et al. | (l2012al) . 

The manuscript is organized as follows: In sect. [IT] we describe the theoretical approach starting from the microscopic 
level and deducing the schematic model later on used to analyse the data. Section UTT1 compares microscopic theoretical 
calculations for hard disks in two dimensions with fits from the schematic model in order to interpret the model 
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parameters and their trends with external control variables like density and shear rate. Sections II VI and [V] describe 
concisely the experimental system and the used Brownian dynamics simulations. In sect. IVI1 we compare schematic 
model and experimental rheology focusing on the stress overshoot as a novel aspect captured theoretically. Finally, 
sect. IVHI mirrors this considering simulation data. A short conclusion summarizes the main findings. 



II. INTEGRATION THROUGH TRANSIENTS APPROACH 



Starting from Smoluchowski dynamics, mode-coupling theory (MCT) uses a Zwanzig-Mori like projection-operator 
formalism and factorization approximations in order to provide a closed set of equations of motion for the intermediate 
scattering function or so-called transient density correlator <I> q (i) for wave vector q. It encodes the slow structural 
rearrangements in the system. Projections on particle-density fluctuations allow us, using $ q (t), to approximate 
ensemble averages like the macroscopic shear stress cr under imposed flow by integrating through the t ransien t (ITT ) 
dynamics. An extensive discussion of quiescent MCT has been given by iGotze and Siogren| (|l992f); iGotze l ([2009) , 
and its ITT generalization to steadily flowing dispersions has been explicitly describ ed bvlFuchs and Catesl (|2009). 
Short reviews with focus on application to colloid experiments have been given by ICrassous et all ( 2008 ) and by 
ISiebenburger et all (|2012h . 

MCT focuses on the dramatic slowing down of the structural dynamics when approaching a critical packing fraction 
ip c . Below, structural correlations relax to zero with an intrinsic final (or a-) relaxation time r Q that far exceeds the 
dilute limit, Brownian diffusion time. Above, the density correlator stays finite even for long times indicating that an 
ideal glass is described. The tp c is thus called the glass transition density. An equivalent discussion holds for the glass 
transition temperature T c . The idealization of a true glass state often does not hold and (littl e understood) activated 
processes turn the glass into an highly viscous liquid (jCrassous et all , 120081 : ICh en et al. .. 2010). In this work however, 
we ignore these activated contributions. 

MCT-ITT describes structural relaxation under an applied shear rate 7. The effect of shear becomes relevant when 
I/7 approaches r a . For smaller shear rates a Newtonian plateau in the long time viscosity 770 is recovered. MCT 
recovers Maxwell's relation 770 oc r a . When the dressed Peclet (or Weissenberg) number Pe = jT a > I, shear thinning 
due to a faster, shear-induced relaxation occurs. In the glass, shear melts the arrested correlators, which makes them 
decay to zero. For 7 — > no Newtonian viscosity can be observed. Instead, a dynamic yield stress tr + = 17(7 —> 0) > 
occurs, which indicates that the Newtonian viscosity is infinite. MCT-ITT predictions should thus be tested in a 
region where the Peclet number Peo = 6irrj s Rflj / (fcsT) <C 1 (rj s solvent viscosity, Rh hydrodynamic particle radius) 
and the Weissenberg number Pe > 1. The former condition is required as MCT-ITT uses the equilibrium structure 
factor S q or more precisely the direct correlation function c q connected to it via the Ornstein-Ze rnicke relation, 
S q = 1/(1 — nc q ) with n the density, to approximate the particle interactions ( Fuchs and Catesl . l2009t ) . 



A. Microscopic MCT-ITT 



The microscopic formulation provides the basis for the understanding of nonlinear rheological properties in MCT- 
ITT. It consists of an equation of motion for the transient density correlator containing as crucial term a retarded 
friction kernel. In the familiar mode coupling approximation, it is expressed as quadratic form in the density correlators 
so that closed equations are obtained. The physical picture behind the approximation is that the fluctuating stresses 
entering the friction kernel participate in the slow structural relaxation and c a n be approxim ated with the lowest order, 
non- vanishing overlap with density fluctuations ( Gotze and Siogreni Il992b iGotzd . 12009). The equations of motion 
are completely specified by the equilibrium structure factor except for one time scale that routinely is matched to 
the short time dynamics. Under shear, the affine deformation decorrelates the stress fluctuations so that the friction 
kernel is forced to decay to zero at long times. Technically, the affine motion enters via an advected wavevector, 
q(i) = (q x ,q y — jtq x ,q z ), which shifts with time and thus tests the effe ctive interactions (i.e. cu ^) a t (for long times) 
increasing wavevector corresponding to shorter and shorter distances ( Fuchs and Catesl . l2002l [20091 ). 

The numerical solution of the microscopic MCT-ITT equations in two spatial dimensions for the correlator, defined 
as the overlap of a density fluctuation at time t with o ne at time t = (wh ose wave vector w a s adv ected in time) 
*q(*) = ($Pq ^Pq(t)(*)) /NS q were presented in detail bv lHenrich et all (|2009l ) and lKriiger et all (|201ll) . The average 
is performed with the equilibrium distribution function so that 3> q (i) is a transient correlator which describes the 
decay of thermally-excited density fluctuations under the action of internal dynamics and external shear. iKriiger et all 
(2011) focused on the incoherent correlator <&*(£) as well as on the mean squared displacements. Even though the 
theory suggests that the correlators in two and three dimensions should be quite similar, and thus the $ q (t) would 
be of interest to experiments, too, we will not repeat the discussion, as the present focus is on developing a simplified 
'schematic' model. The transient correlators are used to find inter alia the shear modulus from an approximate, 



4 



microscopic expression containing information about the shear-driven and internal structural relaxation via $ q (i)- 
It is this relation that we focus on in order to develop a more realistic schematic model, and thus we recall it 
( Fuchs and Catesl . l2009f ). Stress fluctuations at time t after switch-on of the shear rate 7 result from distortions of 



the structure encoded in the (transient) shear-distortion of the structure factor 5S^(t; 7) = Sk(t;7) — Sk- 

^W = ^/t^7^4W;7), (i) 



(27T)° 

where the direct correlation function plays the role of an effective potential (with c' k = dct/dk and particle density 
n = N/V). While the relation is in Fourier-space, as are all microscopic MCT relations (to profit from translational 
invariance) , it bears clos e similarity to approache s that calculate the shear-distorted pair correlation function, and de- 



rive the stresses from it ( Brady and Morrisl . ll997l) . The appearance of the direct correlation function, which in density 



functional theory describes the free energy increase quadratic in density fluc tuations, suggests that M CT implicitly 
assum es a Gaussian distribution of the time-dependent density fluctuations ( Fuchs and Catesl . 120091 : iBrader et all 



2008). The distorted structure factor arises from the time integral over the affinely advected equilibrium structure 
factor weighted by the time dependence of the transient density correlator. The latter captures the memory stored in 
the system of the distortion of Sk at the earlier times. 

pt 

ssu(t-n) = jf di'_g^>$£ ( _ t , )( i'). (2) 

This equation suggests a decomposition of the structural change into a plastic deformation described by the decay of the 
correlator and a (nonlinear) affine distortion described by the wave vector advection in Sk(~t) ■ Quit e intuitively, if no 
struct ural relaxation happens, which would correspond to an purely elastic material, $q(t) = 1 that I Voigtmann et all 
(120121) call 'anelastic', the affine distortion wipes away structural correlations at long times, S^ ncl (t — > 00; 7) — > 1. 



This structural change is elastic, as reversing the flow at some intermediate time, would reverse the effect. In the 
general case, structural relaxation leads to either plastic deformation, a solid-like response when some frozen-in 
structural component remains, $ q (t — > 00) > 0, or to plastic flow, when the structural correlations decay to zero, 
&q(t — > 00) = 0. In consequence, eq. ([I} leads to either anelastic or plastic stress decay arising from the competition 
of structural memory encoded in the $ q (i) and affine motion. 

In startup flows, the above results can be rewritten in the form of an approximated (generalized) Green-Kubo 
relation 



t 



<r(t,7)=7 / dt'g(t',j), (3) 
Jo 

where a shear- and time-dependent stress correlator 5(^,7) captures the nonlinear response. Equations (fTJ) and @ 
give for it: 



k B T f dk 



klky(-t)ky 'Sfc'Sfc(-t) 

kk(-t) Si 



*£(-*)(*)■ (4) 



When its shear-depen dence can be neglected, the familiar MCT approximation for the quiescent stress autocorrelation 
function is recovered (|Gotzd.l2009h . and for finite shear it is a functional quadratic in the density correlator with so- 
called vertices that depend on time via the advected wave vector. The flow curve, viz. the stationary relation 17(7) is 
obtained from eq. [3]) in the long time limit, a(j) = a(t 00, 7). 

Equations ((3]) and (|4|) form the basis for our development of a new schematic model in the next section. Afterwards 
we compare it with explicit numerical solutions of eqs. ([3]) and (f5]) for hard disks in two dimensions in order to gain 
more insight into the parameters of the schematic model from the microscopic approach. Finally, it is tested with 
data from experiment and simulations. 



B. Schematic MCT for viscoelastic media 



In this section we develop and describe a schematic model neglecting the wavevector dependence, i.e. $ q (t) — > 3>(i), 
while keeping the form and asymptotic behaviour near the glass transition of the equations of motion of < f , q (t). Our 
aim is to explain flow curves, frequency dependent linear stress response moduli, and (start-up) stress as function of 
(accumulated) strain for given shear rate switched on at time t = 0. The focus here is on t h e latter, as schematic MCT 
for the former two has been explained and applied to experiments bv lCrassous et al. I (|2008h : ISiebenbiireer" et al. (2009). 
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To br iefly summarize, the equation of motion of the density correlator used previously reads ( Fuchs and Catesl 
2003) 



dt$(t) + T j$(i) + J dt'm(t - t')d' t $(t') 
v 1 $>(t) + v 2 <$> 2 (t) 



= 0, with 



m(t) 



i + {it/icf 



(5) 
(6) 



Short time Brownian motion initiates $(t) = 1 — Yt. The initial decay rate Y describes microscopic dynamics on 
short time scales and depends on structural and hydr odynamic correlations, alb eit MCT neglects the latter. As a 
technical aside let us note already here that, following ISiebenbiirger et all ( 20091 ). we will report times measured in 
the Brownian units D^/R 2 H of the experiments, see sect. IIV1 and not in order to facilitate comparison with 

experiments. The memory kernel m(t) captures the dynamic arrest of slow density fluctuations when approaching the 
volume fraction of the glass transition from below. It prevents &( t) from decay ing to zero in the glassy phase (ideal 
glass). The memory kernel of the repeatedly used F\2 model by iGotzel (1984) is taken. An additional strain term 
jt/j c in the denominator of m(t), with constant shear rate 7, accounts for the shear induced decay of the transient 

correlator to zero (s hear melting) . This shear-rate generalization was called the model bv lFuchs and Catesl ( 20031 ) 
and is reviewed by iFuchsl ( 2010() . The model's parameter 7 C sets the influence shear has on the memory kernel, or 
more precisely, the strain scale for shear-driven correlator decay. The two so-called vertices v\ — v\ + eA/(l — A), 
with v\ = (2 A — 1)/A 2 , and v 2 — v 2 — 1/A 2 , with 1/2 ^ A ^ 1, describe a bifurcation scenario of the non-ergodicity 
parameter f(e) = $(t -t 00, e), which obeys f(e < 0) = and f(e > 0,7 = 0) > /(0,7 = 0) > (type-B glass 
transition). The separation parameter e gives the distance to the transition. We also call / the glass form factor, 
as it describes the frozen-in structure of the glass, which only exists in the quiescent state, 7 = 0. The quantity 
f c = f(e — 0) is additionally called the plateau value, as it gives the amplitude of the final a-relaxation in the fluid 
(asymptotically) close to the bifurcation. Often A = l/y/2, i.e. f c = (1 — A) ~ 0.293, is chosen, which we adopt as well. 
This choice is the same as in prev i ous ap pl ications of this mod e l, whe r e experimenta ll y obt ained stress measurements 
were described bv lCrassous et all ( 2008h : ISiebenbiirger et oil (|2009h : iBrader et all ( 2010h . A is chosen to make the 
asymptotic be haviour (e — > 0~) o f the a - decay time T n o f similar to the one obtained from microscopic MCT for 
hard spheres (jGotze and Sjogren! Il992t ICrassous et all , [200*1. which yields r a oc (— e)" 7 , with 7 = 2.46. Here we 
choose 7(A = l/y/2) = 2.34. 

The numerically calculated correlator $(£) is used in our nonlinear constitutive equation to describe the shear stress 
a(t) (tensorial form neglected) of a colloidal dispersion 



t rt 

dt'g(t-t',[y])i(t')=i dsg{s,j), with 



7 (t, 7) = Va(t, 7) $ 2 (^, 7) + Voc S(t - 0+). 



(7) 
(8) 



In our case, constant shearing starts at t = 0. The generalized shear modulus g(t,j) is chosen as a quadratic form 
in the correlator to mimic eq. (j4])- A function V a (t, 7) characterizes the strength of fluctuations and corresponds to 
the microscopic vertex made up of the static structure factor evaluated at the advected wave vector (viz. the square 
bracket in eq.Q). The time- and shear-rate dependence of ^(t, 7) is the central generalization in our work going 



beyond previous applications of the model by e.g. ICrassous et Ml (|2008l ); ISiebenbiirger et all (|2009l ). where a 



constant v* = v a (t = 0,7) was chosen instead. A short time, high-frequency viscosity, 7700, is added in an ad- hoc 
way in order to account for viscou s proc esses that require no structural relaxation, like the viscosity of the solvent or 
hydrodynamic interactions (jFuchsl . l2010t) . 

The frequency-dependent storage and loss moduli of linear response, G' and G" respectively, are calculated via 
Fourier transforming, while setting 7 = in eq. (|6|) and ([8]). They read 



G?(u)+iG"(u) = iuj / dte- iuJt g(t,j = 0) 



dte^vl $ 2 (<)| 



(9) 



It follows that the plateau value / determines the elastic shear modulus (viz. elastic constant) Goo = G'(lu — > 0,e ^ 
0) = v* f 2 of a glass state. For the fluid, this yields an elastic shear modulus = G'(l/r Q <C oj -C Y, e = 0) as long 
as uJT a ^> 1. 

With constant Vn( t-, A i) = vZ one can f i t flow curves, i.e. 17(7) = a(t — > 00,7) and linear response moduli as shown 



most convinc ingly bv[Siebenburger et all ( 20091 ). Also nonlinear response to large-amplitude oscillatory strain could be 
described by IBrader et all ( 2010|). but will not be a ddressed here. The capability of microscopic MCT- ITT to explain 



stress overshoots was shown by Zausch et all ( 20081 ) within an isotropically sheared hard sphere approximation which 
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however remains numerically demanding. However, the schematic model with constant v* does not capture the stress 
overshoot. Development and validation of a handy schematic model which can recover the richness of stress-strain 
curves in startup flows appears useful and is our main topic. 

We therefore generalize the model to mimic experimentally achieved shear-stress vs strain curves, while con- 
sistently keeping its features and fit parameters used for flow curves and linear moduli. To have a maximum of cr(t), 
i.e. a stress overshoot, one needs the generalized shear modulus g(t,j) to take negative values, cf. eq. j7]). Only 
v,jit, 7) can change sign, since the quadratic form $(i) 2 is inherent in MCT, cf. eq. For enough k values, the 
product of structure-factor deviations S' k S' k ^_ t ^ becomes negative when t approaches r Q causing the A;-space integral 
to become negative. This has been discussed in detail bv lZausch et all (2008). In the sche matic model, Vnit.i) takes 
this role and must be chosen properly. It should be simple, motivated by the results bv IZausch et all (2.008), and 



keep all (desirable) features of the f[?) model, like having a certain set of model parameters for a certain set of 
experimental parameters (temperature, etc.). Note at this point that the present work is not a continuation of the 



stress-tensorial generalizations done to the F$ model in iBrader et Ml ()2009f) . but rather a generalization regarding 



the time dependence of the vertex. 



* b 



b 




log w (W /R 



Figure 1. a): The tim e-dependent vertex v a (t, 7)/^* , eq. ||TDJ, is plotted to show that it mimics the results from the approximate 
microscopic model bv lZausch et all l|2008l ): it agrees qualitatively with fig. 11 there. The model parameters are taken from a 
comparison to experiments (tab.[ll]l8 C (fluid phase)), and from left to right the shear rate decreases, Peo = 4.03- io - ^ 5,6,7,8 ^ . 
b): The solid curves give the corresponding rescaled generalized shear modulus <?(f,7) for the same model parameters. The 
black bars mark the elastic constants G^/vJ = f% (lower bar) and Goc/v% = f(e = 1CP 4 ) 2 (upper bar). The inset shows g(t,7) 
vs strain. The lowest curve is close to equilibrium and the shear-modulus undershoot is about to vanish, cf. the inset of fig. [3] 
b), where this undershoot is completely absent. The dash-dotted curves are glass curves with the parameters from tab. HT1 15°C 
(glass phase), and Peo = 4.03 • 10~ 8 . 



We suggest as form of u CT (i,7) 



v a (t,j) = < (l - (7V7*) 4 ) exp [- (t*/7~) 



(10) 



The strain parameter 7* gives the zero of the stress vertex v a (t, 7), which translates into the zero of the shear modulus 
g(t,j). It thus locates the peak of the stress overshoot and can be read off directly from experimental stress-strain 
data. The strain parameter 7** enforces the final decay of the vertex and thus sets an upper limit to the decay of the 
stress-strain curve to the long time value, which gives the ordinates of the flow curve. These two model parameters 
are regarded as m aterial constants connected to the melting/cage-breaking process, cf. the Lindemann criterion 
(jLindemannl . Il910l ). 

For fixed 7* and 7** and in the asymptotic regime Peo <C 1 and |e| — > 0, our v a (t, 7) form is in full accordance 
with the strain dependence of g(t, 7) in microscopic MCT (see sect. Ill Al especially fig. 5). If the time-scale separation 
to (diffusive) local motion breaks down in $k(-t)(*)i f ne weights of the purely strain dependent vertices in eq. (|3|) 
change. This results effectively in a shear rate dependence of 7* and 7**. We then have to fit these values to each 
experimental stress overshoot, and require the parameters to approach limiting values for decreasing shear rate. 



7 



As implemented below, the choice of eq. (ITU)) for v„(t, 7) is able to satisfy three important requirements, while 
producing quite reasonable fits, see fig. [5] Those three requirements are first monotonicity of the flow curves (some 
are shown in figs. |4] and [8]). Second, a not too long and large valued interval of negative V a (t, 7) is needed to prevent 
a(t) from becoming negative; this is managed by the exponential. Third, one would prefer not to change fits of flow 
curves and linear moduli when fitting the stress overshoot, while of course preferably few new parameters should be 
introduced. Points one and two are not assured by the form of v a {t, 7) itself but only by the procedure described later 
to match its parameters to data. Actually, the negative portion introduced by fcr(7*/7 < t) < into the generalized 
shear modulus 3(^,7) leads for fixed 7* and 7** to a nonmonotonic flow curve at higher Peclet numbers, outside the 
range of validity of the schematic model simplifications. And, for too strong negative tails in u CT (£,7), (unphysical) 
negative stresses may result, which again indicates t he limitation of the m odel. 

The form of v (T (i,7) is motivated by the results of lZausch et all (2008). To demonstrate this, fig. [TJ shows plots of 
eqs. ([8]) and (fT0|k shear rates are given in Peclet numbers, Peg = R^^/Dq, corresponding to the microgel experiments 
from Sect. V in order to ease the latter comparison. The shapes of the moduli undershoots are in good agreement 
with fig. 11 of lZausch et aL (|2008). Al s o, the effect v a (t,j) has on g(t,j) is illustrated, cf. the features of microscopic 
MCT drawn in fig. 11 of lZausch et all 0008). Here and in Sects. VI and VII, the parameters 7* and 7** of eq. (fT0[) 
will be handled as fitting parameters for individual strain-stress curves of different shear rates in order to investigate 
these new parameters in detail. Besides these 7-dependent parameters, shear rates are compared under the same 
experimental conditions (density, temperature, etc.) and shall therefore be modeled with the same parameters v*, 
r, e, 7 C and rj^. The polynomial term in eq. (fTUf describes fast decay of 5(^,7) below zero at 7*. This models 
well the form of a stress overshoot as function of time. The curvature of the final relaxation on the stress plateau is 
modeled by the compressed exponential, as a simple exponential decays too slowly when compared to experiments. 
The parameter 7** serves to adjust the final stress plateau as explained below. The curve of lowest shear rate in fig. [TJ 
has a (nearly) vanishing g(t, 7)-undershoot leading to a (nearly) vanishing stress overshoot, because 7*/7 exceeds the 
t q and &(t) decays to zero before v cr (i,7) causes a stress overshoot. So theory predicts a vanishing stress overshoot 
in the fluid phase for low enough shear rates, which agrees with the constraints of linear response th eory, where i t 
can be shown that the equilibrium shear modulus g cq (t) — g(t,j = 0,e < 0) is completely monotone (Gotze, 2009). 
As expected, this is the c ase in our experiments when Pe becomes smaller t han unity, see fig. [5] a) , while simulations 
bv iFurukawa et all ( 20091 ) and a nonlinear Langevin equation approach bv ISaltzman et aZ.I (|2008h found it only for 
smaller Pe values still. 

The exponents 4 of the polynomial and compressed exponential in eq. (JTUJ are chosen to produce the best fits 
to experimentally measured stress-strain curves. It is necessary to choose even exponents, as strain rate inversion 
should produce the same stress-curve form. A comparison of models with different exponents is shown in fig. O What 
one needs is a function that decays fast below zero near the shear melting point, i.e. (7^/7*)° and a function that 
pulls it back to zero fast enough to prevent a{t) from decaying too much, i.e. exp [—(7^/7**)''] • A choice a = b = 2 
would yield the same qualitative results, but quantitatively fit not as well the experimental data. Higher b exponents 
accelerate the final decay of a to the flow curve value. Higher values of a steepens the slope of Vo-(t,7) around zero. 
Raising both exponents therefore has the effect of amplifying the stress overshoot. 



III. VALIDATION OF THE SCHEMATIC MODEL IN THE MICROSCOPIC MCT-ITT 



Equations Q and (UJ) of the microscopic MCT-ITT were the starting point for the development of our schematic 
model with a time-dependent stress vertex. Its fits to experimental data in sect. I VII will reveal a strong coupling of 
all strain parameters that entered the schematic model, cf. fig. 1101 In this section we compare the schematic model 
with explicit numerical solutions of the wave vector dependent eqs. ([3]) and (UJ) for hard disks in two dimensions. By 
this we will establish that the relations between the model parameters are supported by the microscopic approach. 
Asymptotically close to the MCT bifurcation, all parameters can be taken as constant except for the separation 
parameter e, which measures the distance to the transition. 

Here we show th e time-dependent moduli and stresses resulting from the microscopic MCT-ITT calculations by 
iKriiger et al. 

I ()201lD . which were not shown by these authors, demonstrating that these, as well as our schematic 
model, capture the phenomena found i n the experim e nts a nd si mulations. As the so lution for a three-dimensional 
system is still numerically too costly, iHenrich et all ( 20091 ) and IKriiger et all ( 201ll ) restricted their work to two 
dimensions, considering a system of equal sized hard disks of diameter d, which sets the unit of length. The only 
thermodynamic con trol parameter is the area fraction tp = ^S-- The spherical grid where $ q (i) was computed by 
Krii ger et all uses a discretization of 100 points in radial directi on, q = 0.2, . 6, 1.0, 38.8. The angular 

space was divided in 96 portions, i.e. 6 q = 0.065, 0.13, . . . , 2n (note that lHenrich~ aH (120091) used a cartesian grid; 
the resulting solutions are very similar). From this discretization follows the critical packing of ip c = 0.6985658 and 
the exponent parameter A = 0.7155. The latter determines all power-law exponents of the theory. These values differ 
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Figure 2. An experimental stress vs strain curve of a glass state from fig. [8] at Peo = 4.03 ■ 10 -4 (green squares) is compared 
with fits using the model for vertices v a (t,j) with different exponents a and b (black curves). For all curves, 7* = 0.129 
was fitted. The dot-dashed curve has vertex parameters a = b — 2, -y c — 1.1, and 7** = 0.154 , the dotted one a — 2, b — 4, 
7c = 0.95, and 7** = 0.200, the dashed one a — 4, b = 2, 7 C = 0.95, and 7** = 0.128, and the solid one a — b — 4, 7 C = 0.65, 
and 7** = 0.180. The other parameters are from tab.[TT]T = 15°C. The inset shows the corresponding generalized shear moduli 
<?(i,7) around the a decay time. 



negligibly from the ones found in the equilibrium study of this system by iBaver et al. I <|2007l) (<p c 
A = 0.7167), who used a finer discretization of g-space. 



0.696810890 and 



A. Comparison of schematic and microscopic MCT-ITT 



Figures [3J and |4] show results of the 2d microscopic MCT calculation described in section III Al The generalized 
time-dependent shear moduli shown in fig. [3J contain the information equivalent to the linear response moduli and 
the stress vs strain curves. The microscopic g(t, 7) are fitted by the schematic model together with flow curves and 
stress-strain curves. First, e is tuned to describe the plateau length of the fluid a-decay in fig. [3] a) (purple curve, 
Peo — 10~ 6 ). Then T is used to fit the timescale and together with v* the (3 decay onto the plateau and its height. 
One can use the flow curve, see fig. S]a), to fit j c preliminarily. Then, 7* is determined from the zeros of 7). 
Asymptotically, for vanishing shear rate and separation parameter e, microscopic MCT finds that the stress-overshoot 
peak positions are the same, cf. fig. H]b)-d) and fig. [5j Finally, 7 C and 7** were fitted iteratively to stress vs strain 
and flow curves. We find that asymptotically near the critical packing fraction, viz. small \ (tp— tp c )/ip c \, all schematic 
model parameters can be taken as constant for the comparison with the microscopic calculations. We thus drop 
the dependence of ^(t, 7) on 7 in this section, because v a (^) depends effectively just on strain 7 = jt. Only, the 
separation parameter e varies linearly with (ip — tp c )/tp c , as required; the chosen parameters of the schematic model 
are given in tab. Q] 

Overall, the schematic model fits the microscopic moduli quite well outside the short time regime. In fig. [3J the 
plateau values are matched well, while only for the large Peo = 10~ 2;-3 , the separation of a- and short-time processes 
deteriorates and systematic deviations around the plateau appear. A closer look at the region where g(t, 7) is negative, 
see the inset in fig. [3Ji), however, reveals another, presently more relevant, systematic deviation. The microscopic 
calculations have a less-compressed a-decay, and a broader negative minimum, than the schematic model. Fig. [3] b) 
shows that a time dependent vertex v a (^f) with two exponents a = b = 2 agrees better with the microscopic MCT-ITT. 
It fits the moduli undershoot and thus stress overshoot better than with two exponents of a = b = 4, as in eq. (fTOf . 
This indicates that, compared to the data from experiments and simulations (cf. sects. PVTl and IVlTj) the microscopic 
MCT-ITT transient correlators decay too slowly in the final step; recall fig. where experimental data agree better 
with ^,7(7) for a = b = 4 than for a = b = 2. This deficiency of the microscopic theory had already been observed 
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Figure 3. Panel a), rescaled generalized shear moduli from microscopic (thick, colored) and schematic (thin, black) MCT. Fluid 
curves are purple ((if — ip c )/ip c = — 1CP 3 ), the glass curve ((tp — ip c )/f c — + ICP 3 ) is green, and Peos are given in the picture. 
For the parameters see tab. U (fluid lines solid, glass line dashed). The exponents in v CT (7) are set to a — b = 4. The 
black bars mark G^/v^ = /,? (lower bar) and Goc/v% = f(e = 8 ■ fO -4 ) 2 (upper bar). The inset shows the final decay in g(t,j) 
vs strain. Small differences between microscopic and schematic model become noticeable. The leftmost curve approaches the 
equilibrium relaxation and thus decays prior to the shear- melting strain 7*; it thus leads to no stress overshoot. 
In panel b), the schematic vertex v a (-y) with both exponents set to a — b = 2 is shown. It fits the microscopic MCT moduli 
better. 



by iKriiger et all ( 201ll ). They found that the transient correlators around their plateau value from the microscopic 
MCT-ITT agree with the ones from the simulations up to a rescaling of the shear rate by a factor around 5. Yet, 
during the final decay the actual decay of transient or stress correlators is faster than theoretically predicted. In total, 
this leads to the strong overestimate of the strain 7* of the stress maximum by the microscopic theory; this will be 
discussed next in connection with fig. [SJ In the schematic model, the freedom in the choice of the strain-dependence 
of iv(7), i-e. via 7* and 7**, can be used to provide better fits. 

Figure [?] shows the corresponding flow curves and stress vs strain curves from the two-dimensional microscopic 
MCT calculations and the best fits of the schematic model to it; the parameters are as used in fig. [3] and listed in 
tab. [T] As discussed bv lHenrich et al. (2009), the flow curves from microscopic and schematic MCT-ITT show a yield 
stress and the characteristic change from a fluid to a solid like stress curve as function of shear rate when crossing 
the (ideal) glass transition. Panel a) of fig. 2] extends the flow curves of the schematic model beyond its range of 
applicability to indicate that it predicts a non-monotonic behavior of the stress for too large shear rates; as the inset 
in this panel shows, for high Peclet numbers, the negative tail of the vertex- function v a (^) from eq. (flTJ)) dominates 
the generalized shear modulus and leads, for large 7, to decreasing stresses, which can even become negative. Clearly, 
the schematic model approximation to use a single vertex- function v a (^) with strong negative tail breaks down here; 
the microscopic MCT-ITT continues to predict a monotonic flow curve. The stress vs strain curves show a stress 
overshoot for all densities and shear rates outside the linear response regime. The linear response regime only exists in 
fluid states, where it can be observed as long as the dressed Peclet number Pe = jr a is small compared to unity. For 
small strains, the stress increases linearly with strain and a common elastic constant G oo can be observed, o~ ^ G-ooT- 
Only for bare Peclet numbers approaching unity, the merging of short and long time processes results in a shift of 
the linear portion in the stress-strain curve. The stress overshoot in the microscopic MCT calculation is broader 
than the one in the schematic model as the latter has a more rapidly varying time dependence at long times. The 
overall satisfactory agreement of the moduli g(t, 7) from the schematic and the microscopic calculation translates into 
a comparable agreement of the stress vs strain curves. Yet, the differences between both models during the final 
decay are more prominent in the stress vs strain curves, where clear deviations in the shape of the stress overshoot are 
noticeable. Changing the exponents in the schematic model stress vertex v CT (7) to a = b = 2 significantly improves the 
agreement of the schematic to the microscopic MCT stress-strain curves, as shown in panel d) of fig. 0] This suggests 
that the actual shape of the stress overshoot could provide a sensitive test of particle properties and interactions in 
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Figure 4. Results from the microscopic MCT of hard disks (symbols) for three packing fractions close to the glass transition 
compared to fits using the schematic Ffjj -model (lines) where only e has been changed. In a), the flow curves are shown, in 
b) the stress vs strain curves of the fluid state (((p — tp c )/tp c = — 10~ 3 , purple circles), in c) the critical ones (tp — cp c , cyan 
triangles), and in d), the stress vs strain data of the glass state ((ip — ip c )/ip c — + 10 , green squares). The black, continuous 
curves are F^ fits with parameters in tab. [T] the panel d) includes stress-strain curves obtained from the v a (7)-vertex with 
two exponents a = b = 2 (as in fig. \Sj))) to show the improved agreement (lower set of curves, which are shifted for clarity). 
The inset in panel a) shows the generalized modulus g(t,-y) from the schematic model as function of linear time for two high 
shear rates, Peo = 0.1 and 1 as labeled. The effect of the time-dependence of f CT (7) is studied; it is included {black solid line), 
or neglected, t> CT (7) = v% (blue dashed-dotted line). For high shear rates, the negative tail in the modulus causes the steady 
stress to decrease as noticeable in the main panel a). 



viscoelastic dispersions. In lKoumakis et al\ ( 2012bl ) stress-overshoots with respect to different interaction potentials 
were investigated. 
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Table I. F^ parameters used to fit the microscopic MCT calculations; the units are [v *] = k B T/d 2 and [r] = D /d 2 



Asymptotically (limit of Pcq <^ 1 and e — > 0), microscopic theory predicts that the transient correlators depend on 
the strain only (time and strain rate reduce to one master variable). Thus, the stress maximum lies at asymptotically 
identical strain values 7*, and a strict relation between 7* and 7** holds. The schematic model fits to the microscopic 
MCT-ITT curves were performed, accordingly. Looking in more detail at the microscopic calculations, this can be 
observed in the glass, e > 0, where the 7*-values become independent on shear rate in fig. [5] Conveniently, the 7* 
are taken from the zeros of the microscopic modulus eq. ([J}. In the glass, with ((p — (p c )/(p c = 10~ 3 , all 7*s up to 
Pe — 10~ 3 are essentially the same. For Pe = 10" 2 , when the strain-rate independent short-time decay becomes 
relevant, an increase of 7* sets in. In the critical case, the j3 decay onto the correlator plateau is stretched to a 
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maximal extent, i.e. reaching constant 7*s requires smaller Peo- In the fluid, the stress overshoot vanishes when the 
linear response regime is approached; the inset of fig. [5] shows the relative peak height 0^/(7(7) — 1 from microscopic 
MCT-ITT. For the five lowest shear rates, no stress overshoot can be observed, which makes 7* not well defined. The 
zero in g(t,j) shifts because of the form-change in the stress vs strain curve but not because of a variation of the 
parameter 7* in the model. We thus take 7* as constant asymptotically in the fluid, also, as is suggested by figs. [3] 
and|4j where it was constant throughout. Outside the asymptotic constant regime, 7* grows monotonically with shear 
rate. 

An important difference between the microscopically calculated maximum strain values 7* and the ones obtained 
from the fits to the experimental and simulational data shall already be mentioned even though it requires to compare 
fig. [3] with the later figs. [TU] and Q2J It concerns a (large) quantitative error of microscopic MCT-ITT. It predicts 
the stress overshoots to arise at around a factor of two larger strain values 7* than seen in the experimental data, 
fig. [TU] inset. While in qualitative agreement — in all cases a constant 7* value is approached for decreasing shear 
rate — comparing the numerical values at the glass transition (at s — + ), gives quite similar estimates from the 
d = 3 experiments (7* w 0.07) and the d — 2 simulations (7* 0.04), which, however, differ appreciably from 
the theoretical calculation (7* ps 0.22). As discussed above, this results mainly from the too slow variation of the 
theoretical correlators during the final relaxation step. The schematic model which can be validated from comparing 
it with the microscopic calculations offers a handy description of experimental data as its parameters can be adjusted 
more freely. 

The relative peak height, max[<r(t)]/cr(oo) — 1 = <7pk/c(7) — 1 can also be calculated from microscopic MCT and 
is shown in the inset of fig. [5] for various packing fractions (ip — ip c )/ip c in a range from — 10~ 2 to +10~ 2 . Note that 
the theory applies close to the glass transition at tp c and makes no predictions for dilute or almost closely packed 
dispersions. For all shear rates, the relative peak height increases with packing fraction linearly. We take this as 
a simple statement of a Taylor-expansion like relation, where the prefactor, however, is surprisingly large so that 
changing packing fraction by 2 % in crease the relativ e peak height by a factor around 2. This result of microscopic 
MCT-ITT contradicts the results by iKoumakis et al\ ( 2012a[ ). where the opposite was measured. A decrease of the 
relative peak height with packing fraction was attri buted to approach ing random close packing. Considering the 
known sensitivity of the stress overshoot on ageing (|Derec et aZ.1 . 120031). the lack of i nforma tion on the preshearing 
protocol let us speculate that too short waiting times were used by Koumakis et al\ ( 2012alh our results correspond 
to the limit of starting from an equilibrated sample as was checked in the sim ulations and could experimentally be 
achieved because of a hopping process melting the glass ( Crassous et all |2008| ). 
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Figure 5. The strain values 7* of the stress maximum from microscopic MCT are given as function of shear rate (rescaled as 
Peo = 'yd 2 /Do) corresponding to the fluid (circles), critical (triangles), and glass (diamonds) states of fig. [4] 7* are read off 
from the zeros of g(t,y), eq. Q. The inset shows the relative peak height a p k/o- 3t — 1 of the stress overshoot from microscopic 
MCT for several Peo and packing fractions tp around the critical one tp c . Peak-heights increase with shear rate and with packing 
fraction. 
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IV. EXPERIMENTAL 



Experiments on the start -up stress strain curv es were performed on a we ll studied mod e l glas s former, whose 
rh eology was investiga t ed by ICrassous et al. | 1)20081) : ISiebenbiirger" et al. | 1)20091) iBrader et al. (2010) and is reviewed 
bv ISiebenburger et al\ (|2012t >. 

For the rheological measurements, we used an aqueous suspension of core-shell latices with a polydispersity of 
17 %, which conveniently suppresses crystallization. The core consists of poly(styrene) (PS) onto which a ther- 
mosensitive crosslinked network of poly(N-isopropylacrylamide) (PNiPAM) is affixed. The degree of the cross-linker 
N,N'-methylenebisacrylamide is 2.5mol%. Below the lower critical solution temperature (LCST), the PNiPAM-shell 
is swollen and shows a line ar shrin king up to 2 5 °C wi th increasing temperature. The same latex was already used 
in ISiebenbiirger et al. (l2009h and IBrader et al\ ( 2010h . and the full characterization of these particles is given in 
ISiebenburger et al ( 20091 ) . Residual charges due to the synthesis were screened with 0.05M KC1, leading to almost 
pure steric interaction. The solid content of the suspension, i.e. the non-evaporable weight percent of the particles 
in the aqueous suspension, was found to be 8.37 wt% ± 0.02 wt%. The effective volume fraction ip e g was calculated 
by using the correlatio n of mass concentra t ion c. hydrodynamic radius Rh and effective volume fraction as shown 
in the inset of fig. 6 in ISiebenbiirger et al. (120091). We expect th e glass transition to lie at ip 9 eS = 0.640 in this (for 
colloids) slightly polydisperse system, cf. Crassous et all ( 20081 ); ISiebenbiirger et ail ( 20091 ). At the temperature of 
15°C a volume fraction <p c ff = 0.65 and for 18°C a (p c g = 0.60 was found. This equates to relative packing fractions 
(fcS — VcffVVeff °f 1-6 ' 10~ 2 and —6.25 • 10~ 2 , respectively. For comparison, the random close packing (RCP) of 
monodisperse hard spheres in 3D is ip rC p = 0.63, but in polydisperse mixtures RCP is known to be much higher, cf. 
ISchaertl and Sillescul (| 19941) . 

The rheological measurements were pe rformed at the rheometer MCR 301 from Anton Paar with a cone-plate 
geometry (further details are described bv ISiebenbiirger et al\ ( 20091 )) with an layer of very low viscous paraffin oil to 
prevent evaporation. All start-up measurements were performed after a pre-shear at a shear rate of 100/s for 200s 
and a waiting time of 10mm. After this waiting time, ageing effects on the stress overshoot are still observable but 
can be considered small in our system even at the glassy packing fraction. Fig. [5] shows the development of the stress 
overshoot for different waiting times up to and beyond the one (t w = Wmin) used later on throughout this work. 
Between lmin and 10mm, the form of the stress overshoot still changes noticeably, while increasing the waiting time 
further to t w = 60min achieves the long-time shape. Then the limit of infinite waiting time assumed in the theory 
is established. In order to sample a large data set, we neglect the small difference between the stress-strain curves at 
t w = Wmin and the completely converged ones; Fig. |5] shows the worst-case di fference in our study. The observed 
magnitude of t w agrees with the results of an aging study in the microg els by C rassous et al. (2008), who argued 
that an ergodicity-restoring ('hopping') process cuts-off aging in the microgel glasses. The frequency test from 10 - 
5T0~ 4 /s was started after a pre-shear at a shear rate of 100/s for 200s and a waiting time of 20s with a logarithmic 
time ramp of 20 - 2000s. The flow curves were measured from 5T0 _5 /s - 1 ■ 10 3 /s with a logarithmic ramp from 
2000s - 20s and back. The bare Peclet numbers were calculated from the viscosity rj s of water and the hydrodynamic 
particle radius according to Peo = 6irr] s R^ 1 'j / (fcgT) . 



V. SIMULATION 



Simulations of start-up stress strain curves were p erformed using a well studi ed two dimensional binary mixture of 
ha rd disks, whose glass tr ansition was d e termi ned bv lWevsser and Ha i nail ( 201ll ) and whose rheology was investigated 
bv lHenrich et al. I (|2009l) : iKriiger et al\ (|201ll) . 

The concept of the algorithm was detaile d by Scala et al. J2 007) in the three dimensional variant. Here we use 
the two dimensional algorithm described bv iHenrich et al\ ( 20091 ) . We consider a binary mixture of hard disks of 
mass mo with the size ratio of S = d s /db — 5/7, where d s and db denote the diameters of the small and big particles 
respectively. N = 1000 hard disks (N s = 500 small and Nb = 500 big ones) move in a two dimensional simulation 
box of volume V at a packing fraction ip = -^(N s d 2 s + A^d 2 ). After placing the particles on their initial positions we 
provide Gaussian distributed velocities with variance ) = v ®i where a denotes the species and i the number of 

the particle. To propagate the system at time t forward in time, we employ a semi-event-driven algorithm. For every 
particle, e.g. for particle i of species a (big or small) at the time t, the algorithm determines the possible collision 
time At"P with any other particle j of species f3. The smallest solution for all particle pairs determines the next event 
in the algorithm. All particles can then be propagated with constant velocity, until at the collision time elastic binary 
collision laws impose new velocities. 

So far, the algorithm described above yields ballistic motion. In order to mimic Brownian motion we modify it by 
introducing a thermostat which at every integer-multiple of the time TBVo/d s = 0.01 triggers a so-called Brownian 
step. In the Brownian step, all particle velocities are freshly drawn from a Gaussian distribution with variance 
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Figure 6. Experimental stress vs strain measurements for different waiting times t w as labeled after a preshearing for 200s 
with Peo = 4.03 • 10 _1 . The Peo for all curves is the same and given in the plot. The curves correspond to the glassy packing 
fraction ip ~ 0.65 adjusted by the same experimental settings as the 15° C curves of fig. H]b). Throughout the work reported 
subsequently, the curves with t w = Wmin are used. 



m o v o IQ^bT) = 2 for all particles. This assures that the particles move diffusively with a short-time diffusion coefficient 
Do/(vod s ) = 0.005 on time scales which are large compared to rg. The Gaussian is centered around zero at each step 
which enforces that the center of mass of the N particles remains fixed. The shear is simply added as y-dependent 
bias in x direction, reading jyt, where 7 defines the shear rate and y\ the y-component of the position of particle i. 
We choose Lees Edwards boundary conditions to simulate a bulk system. 

The system starts from a cubic lattice with all particles having the same size at low density. After a short initial 
equilibration the particles are swollen to the desired packing fraction guaranteeing an amorphous system. It is 
necessary to wait for the system to relax before meaningful averages from the equilibrium distribution can be taken. 
750 independent systems were created for a li quid state at (p = 0.79 and f or a glassy state at <p = 0.81. As we expect 
the glass transition to lie at <Ps im ~ 0.795, cf. IWevsser and Hainall ( 201ll ). this equates to relative packing fractions 



(ip — ^sinJ/Vsim °f — 6.3 ■ 10 3 and 1.9 • 10 2 , respectively. For compa rison, RCP of mono disperse hard disks in 2D is 
frcp = 0.82 and for this binary mixture it is approximately 0.817, cf. lHainal et al\ ( 2009h . 



Equilibration was performed with Newtonian dynamics (without imposing the Brownian step) for 10 5 time steps 
(in units of d s /vo) in the liquid state. In the glassy state 2 • 10 6 time steps were used for equilibration. We assume 
that the system is equilibrated, when the time-depend ent correlation functi o ns do not depend on the time origin. For 
a more detailed explanation for the glassy system see lWevsser and Hainall ( 201ll) . Using these equilibrated systems 
we switch on the shear and generate the transient correlation functions. 

Rece ntly, transient shear b anding has been discussed as being closely connected to an overshoot in the start-up 
stress ( Moorcroft et all 1201 lh ■ We therefore recorded the velocity profiles in the simulation for various shear rates 



and strains; see fig. [7J As the linear velocity profile is added to the gaussian velocities as a bias at every Brownian 
timestep, measurements of the velocities were taken shortly before the Brownian timestep is applied, thus giving the 
maximum time possible to a potentially developing shearband. In total we used 600 independent configurations at 
ip = 0.81 in the glass to determine the velocities in the linear response, stress overshoot and stationary regime. The 
scatter of the measurements around the nominal shear rate gives no indication of shear-banding in the simulation. 
The bias added after the Brownian timestep means that we can rule out slowly developing transient shear bands as 
phenomenon occurring together with stress overshoots in our simulations. 
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Figure 7. Velocity profiles for a) Peo — 0.02 and b) Peo = 0.002 at tp = 0.81 in the glass. The profiles were measured according 
to the procedure described in the text. The different strain rates in the two panels correspond to the linear response regime 
(with red crosses and strain jt = 0.02), the stress overshoot (with green circles and strain jt = 0.08, yt = 0.05) and the 
stationary regime (with blue stars and strain yt = 1.0). The lines with colors, corresponding to the ones of the data points, are 
fits of a linear function to the data. 



VI. COMPARISON OF EXPERIMENT AND MCT 



In the following our new time-dependent v a (t,j) is tested with model experimental data me asured on the ther- 
mose nsitive core-shell particles immersed in water at 18°C (fluid state) and 15°C (glassy state) (jSiebenburger et~aj\ , 
2009). Because the time-dependence of the vertex v a (t, 7) includes negative portions, care has to be taken in matching 
the new schematic model to experimental data. As it is naturally demanding to go deeply into the asymptotic regime 
of MCT's glass transition (i.e. Peo <§C 1 and e — > 0), it is convenient to effectively account for short-time motion 
effects by a explicit shear-rate dependence of the vertex (in addition to a pure strain dependence), which means to 
match 7* and 7** to every shear rate, cf. the discussions in sects. Ill Bl and Mil We found the following procedure 
useful. 



A. Matching the schematic model to experimental data 



The parameters e, 77^, v*, T, and y c , together with a constant stress vertex, v„{t , 7) = vZ, are used to fit the 
flow curve (7(7) and linear moduli G'(u>) and G"(u>), as discussed bv lHajnal and Fuch s (2009). In rescaling time and 
shear rate according to t t— > Tt and 7 t— > j/(Ty c ), one can verify that e and 7700 set the shapes of flow and moduli 
curves, so they are adjusted in a first step. T and v* shift those curves horizontally and vertically, while y c shifts 
only the flow curve diagonally, so these three are adjusted in a second step. One may iterate this to improve fit 
quality. Whenever there is a plateau in G', this can be used to determine v*, because eq. © yields G' = v%f 2 for 
sufficiently long plateaus in the correlator $(t); this requires small Peo (lin. response 7 — > 0) and r a large compared 
to 1/uj so that the fluid reacts elastically. The results are shown in fig. [5] While flow curves and linear- response 
storage modulus G' are described well, the loss modulus G" exhibits deviations at low frequencies. This is ca used by 
a low -frequency process, often called the hopping process, cf. fig-Eflb), w hich is not part of (idealized) MCT ( Gotzd . 
2009); its influence on the spectra was also seen bv lCrassous et al\ { 20081 ). The theory describes an ideal glass, whose 
quiescent correlators $ q (t) do not relax to zero. The measurements indicate that there is a very slow process in the 
glass which leads to a finite Newtonian viscosity far outside of the range accessible by experiment. As an outlook to 
the stress-strain curves let us add that as consequence of this hopping process the measured stress-strain curve at 
7 = 10~ 4 /s in fig. [9]b) possesses no stress overshoot, although MCT predicts it. MCT considers this state to be a 
glass, while the stress strain curve approaches the linear response one characteristic of a fluid. The disagreement is 
again caused by the neglected hopping process. 

From here on, e, rjoo, v% and T are not changed, and the acquired linear moduli are kept. Yet, y c is increased, 
increasing the positive area under the shear modulus g{t, 7), because turning on the time-dependence of V(r(t,j) 
introduces a negative contribution. The aim of the procedure is to leave unchanged the acquired flow curves when, 
now, addressing the start-up stress using the time-dependence of v a (t, 7). For each curve of a set, the peak position of 
stress vs strain is taken as 7*. Where this is not available, 7* is extrapolated as described in the next section. Now, 
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Figure 8. Panel a) shows flow curves and panel b) linear response moduli. Purple circles are experimental measurements at 
18°C (fluid phase) and green squares at 15°C (glassy phase). The upper moduli curve is G' , the lower one G" . The black, 
continuous curves are fits, see tab. [II] for the fit parameters. In the glass phase, the low frequency part of G" can not be 
fitted by MCT and arises from a ('hopping') process that deviates from ideal-glass behaviour l|Siebenburger et aU [2009). The 
dashed horizontal lines mark the elastic constants, g(t,j) = G%o in the fluid phase and G OCI {e = 1CP 4 ) in the glassy phase; they 
give the plateau values of g(t, 7) prior to a decay or shear induced decay, respectively. 



7 C is increased to an appropriate value, i.e. one that approximates the stress curve around its peak. One value of j c 
applies for the whole shear-rate set. Finally, 7** is used to tune the final approach of stress to its long time value on 
the flow curve. This finishes the fitting, which, as described, works step by step. We found this a useful feature, as it 
does not destroy previous achievements, conserves the monotonicity of the flow curve by means of its implementation, 
and is therefore in accordance with previous applications of the schematic model. 




Figure 9. Stress strain curves: Purple disks show experimental data at 18°C (fluid state) and green squares at 15°C (glassy 
state), while black, continuous curves are fits, with parameters listed in tab. [II] for 7* and 7** see fig. 1101 In a), start-up 
viscosities, viz. rj(t, 7) = u(t, 7)/7, are shown because this splays the curves; in this manner of plotting, elastic regimes do 
not coincide as in b), where <r(t, 7) is shown directly. The dashed lines mark the elastic response, a(t) — G^yt (in a)) and 
Goo(e = 10 -4 ) (in b)), with elastic constants taken from fig. [H]b). 

Fig. [9]shows our results of the newly introduced v a (t, 7) in fitting stress-strain curves and stress overshoots measured 
in the core-shell particle dispersions. Panel a) shows <j(i)/7 corresponding to the start-up viscosity in order to splay 
out the curves. The abscissa is given by the accumulated strain 7 = jt. Panel b) gives a(t) itself in order to indicate 
the overlap of the curves in the initial elastic region. The dashed lines in both panels mark the elastic behavior 
expected from the analysis of the linear response moduli, viz. a(t) — G 00(5)7 = v*f 2 ^. The plateau modulus was 
included as horizontal dashed line in fig. [5] This indicates the agreement of the elastic behavior from frequency 
dependent linear response and the linear regime in the stress strain curve. Importantly, the elastic response is given 
by Goo = g cq (t —> 00), which MCT calculates as characteristic for the glass state, and not by g(t = 0,7). This latter, 
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T = 18° C, tp eft = 0.60 


v* = 85 


r = 105 


: = -2 ■ 10 ' 


r?oo = 0.28 


7c = 0.65 


log ln (Pe /4.03) 

olU \ *-* / / 




-8 


-7 


-6 


-5 


-4 


-3 


-2 


-1 


7* 




0.0505 


0.0593 


0.0762 


0.0913 


0.117 


0.153 


0.210 


0.318 




7" = 1.55 7* - 0.0163 


T = 15°C, ip ef f = 0.65 


«* = 123 


r = 125 


£ = 1 • 10" 4 


??oo = 0.30 


7c = 0.65 


log 10 (Pe /4.03) 


-9 


-8 


-7 


-6 


-5 


-4 


-3 


-2 


-1 


7* 


0.0702 


0.0710 


0.0748 


0.0826 


0.104 


0.129 


0.165 


0.227 


0.341 




7** = 1.59 7* - 0.0245 



Table II. Parameters of the schematic model in the fits to the experiments on the core-shell particle dispersions; the units are 
K] = k B T/R 3 H ,[r] = Do/R 2 H , and [ Voo ] = k B T/(D R H ). 



equal-time value of the stress auto-correlation function characterizes stress fluctuations at such short time scale that 
fluid and glass cannot be discerned. In equilibrium it can be calculated by the Bixon-Zwanzig formula, and it also 
does not change when crossing the glass transition, while Goo is finite only in the glass. 

Beyond the linear regime, the stress strain curves exhibit the overshoot. The schematic model captures the non- 
linearity well in position, height and shape. In the fluid phase, the low 7 stress-strain curves, fig. ||a), show the 
theoretically expected behavior, that the stress overshoot vanishes when approaching linear response; this holds as 
the equilibrium stress correlation function g cq (t) is positive and monotone. A clear decrease of stress overshoot 
with decreasing shear r ate in a fluid PMMA (polymethyl methacrylate) hard sphere sample was also observed in 
iKoumakis et al\ ( 2012bl ). Quantitatively, the experimental overshoot vanishes somewhat faster than theory predicts, 
which we blame on the oversimplification of the schematic model. The approach of the experimental data to the 
linear-response behavior in the glass state arises because of the mentioned hopping process, which MCT neglects; the 
schematic model fits consequently show a weak overshoot remaining asymptotically for vanishing shear rate in the 
glass. Theory and experiment give a weak increase of the stress overshoot with shear rate outside the linear response 
region; see sect. Mil and especially fig. 2] for a more detailed theoretical discussion. Overall, we judge the agreement 
between the schematic model and experimental data as good considering the simplification to condense the k-space 
integration of the microscopic MCT in eq. 2] to three schematic model parameters. 

For small strains, the initial decay of <&(t) matters as described by eq. ((5|). It is strain-independent, and yields 
different offsets for MCT's stress-strain curves in the limit of 7 = jt — > 0. For the higher shear rates in fig. |H1 
experiments cannot resolve the small strain windows any more. How well our theory describes the low strain limit in 
this case cannot be checked due to lack of measurement data, which caused us to cut our curves there. 

In the next section, the fitting parameters 7* and 7** are analyzed. Our conclusion is that they are coupled and 
could effectively be reduced to a linear coupling that holds for the whole shear-rate set. 



B. Analysis of overshoot drift in 7* and 7** 



When analyzing data with the new schematic model, 7* can be read off as the strain value where the stress is 
maximal, while 7** is adjusted to retain the fitting of the flow curves, see eq. (|10[) . The measured 7*-values drift 
slightly with shear rate, i.e. the stress peak position 7* (7) is a weakly increasing function of shear rate. The drift 
in 7* over Peo is plotted in the insets of fig. [TO] The dependence of 7* on Peo is quite weak, and, as we argue in 
sect. Mil is due to preasymptotic corrections for shear rates that are outside the asymptotic regime, Peo C 1. A cubic 
dependence of 7* on the logarithm of Peo can provide a guide to the eye through the fitted 7*-values. An asymptotic 
value 7* = 6.98 • 1CT 2 for log 10 Pe Q < -8.13 is suggested for the glassy data. For the fluid phase, 7* for lower Peo 
can not be determined from the fits, because the typical shear melting timescale I/7 then exceeds r Q of the fluid, i.e. 
7* cannot be chosen uniquely. We consider the two lowest 7* points to be effected by this. 

Fig. 1101 shows a plot of 7** vs 7*. Filled symbols give parameters directly taken from the experiment (viz. 7*), or 
well determined by the fit (viz. 7**). The relation between both fit parameters is clearly linear (plus constant). From 
the trend of the black (7* measured) dots in fig. [9j one concludes that 7** really is coupled to 7* and refers to the 
same material parameter, which however exhibits a weak drift dependent on 7. This effectively reduces our two new 
fitting constants to one material parameter, namely 7*, which can be read off or extrapolated from experiment within 
the schematic approach. The second parameter then follows from the linear relation. In glass and fluid states, the 
relation between 7* and 7** is identical within the scatter of the fitted parameters. As consequence, we adjusted those 
{7*, 7**} points that were not extractable from measurement to be on the same linear trend, indicated by straight 
lines in fig. grey symbols there indicate extrapolated parameters. For those, not 7* but its relation to 7** is set, 
and only 7* is fitted. 
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Figure 10. Shown are the strains 7* of the stress maximum and the strains 7** marking the approach to the stationary flow 
curve used in the fits of fig.[9]a) (disks, i.e. fluid phase) and b) (diamonds, i.e. glass phase). The black symbols refer to measured 
7* which suggest a linear relation between 7* and 7**. Grey symbols are extrapolated with this linear trend. The inset shows 
the position of the stress maximum 7*s as function of Peo, with a guide to the eye (cubic polynomial in log Peo). Tab. [TT] 
summarizes the parameters. 



VII. COMPARISON OF SIMULATION AND MCT 

Fig. [HI and [T2l show flow and stress-strain curves from 2d simulations of Brownian hard disks for given shear rates 
and for two densities. Both are very close (below and above, respectiv ely) to the critical density, ip c — 0.796, where 
MCT locates an idealized glass transition ( Wevsser and Hainall . l201lh . Like in the experiments, however, hopping 
processes also melt the computer glass. Now the Peclet number is taken as Peo = 7 <P s /Dq, see section [V] for a 
detailed description of the simulations including of d s and Dq. Frequency dependent linear response moduli are not 
yet a vailable. Therefore, after fitting e to the flow curves (setting 77^ = 0, as there is no solvent ( Henrich et all 
2009)), v* must be read off in fitting the linear, elastic increase of a(t) for small strains. This makes it necessary to 
iterate between flow and stress-strain curve, until the best fit is reached. 

Overall, the schematic model can well describe the simulation data for the typical liquid and typical (MCT-) glass 
state. For Peo £ 0.02; 0.002 our MCT's flow curve does not perfectly fit the simulated one; see fig. [TTJ Simulation 
in this interval has deviations from theory to positive stress values. These cannot be fitted perfectly with a convex 
function, which our theoretical glass flow curve is on this interval. The scatter of the corresponding stress-strain 
curves from simulation (see fig. 1 12[) prevents us from reaching a conclusion on the significance of this deviation. The 
fits to the stress overshoot again favour the form of eq. (fTU|) with its two exponents of 4. It provides the best fits 
to the stress overshoots (alternative fits are not shown but contain a broader and less well defined maximum in the 
stress.). Fig. [T2l shows the outcome. Because the simulations extend to short times, an apparent variation of the 
elastic regime with shear rate (viz. Peclet number) is noticeable. This effect, which, we expect, exists also in the 
experiments but cannot be resolved there, should not be interpreted as a dependence of the elastic shear modulus 
on shear rate 7. Rather, it results from the lack of a clear separation of time scales; see fig. \T$). For (bare) Peclet 
numbers Peo — jd^ / Do approaching unity, the stress correlation function has not completely decayed onto the elastic 
plateau before the time Thus its short time dynamics gets affected by shear, and the stress varies with shear 

rate before the overshoot in fig. [T^l This rather strong effect in the simulations of hard disks may be specific to the 
excluded volume interaction. A direct comparison to experiment is not possible, because there short-time data for 
such high Peos is lacking, cf. fig [HI 

The observation made when analyzing the linear correlation of 7* and 7** in the experiment, cf. fig. 1101 is also 
strengthened by the simulation. For the latter, 7* and 7** show a strong linear correlation, too, see fig.[I31 The linear 
trend slightly differs for the two different densities, which appears reasonable, as packing fractions not very far from 
random close packing were simulated. Generally, it is expected that random close packing causes anomalous stress 
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Figure 11. Flow curves of a Brownian dynamics simulation of a two-dimensional, binary mixture of hard disks at packing 
fractions ip = 0.79 (fluid state,purp/e disks) and ip = 0.81 (glassy state, green squares) are compared to fits using the F^P model 
(black solid lines); parameters are given in tab. IIIII 




Figure 12. The stress vs strain curves from the Brownian dynamics simulation are shown which correspond to the flow curves 
from fig. 1111 Disks in a) are for ip = 0.79, and squares in b) for ip = 0.81; shear rates, converted to Peclet numbers, are as 
labeled. The black lines are fits using the Ffp model with parameters in tab. IIIII 



fluctuations with strong dependences on the control parameters like density. Reassuringly, the linear relation between 
the two strain parameters of the model seems to be a very stable feature of our F^p model for stress overshoots. 
This highlights that a single strain parameter determines the nonlinear rheology and may be considered a material 
constant. We take it to be 7* which is the strain of the maximal stress, which moreover varies only little as function 
of shear rate. Also the simulations suggest that 7* increases with Peo, see fig. 1131 Inset where cubic polynomials in 
log 10 Peo provide guides to the eye. 



VIII. CONCLUSIONS 

We generalized a schematic model of the mode coupling theory of flowing glass formers, the so-called F^p model, 
such that consistently flow curves, linear response shear moduli and stress-strain curves (with a focus on the stress 
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Figure 13. Plot of the strain parameters 7*s from the fits to the simulated stress-strain curves in fig. 1121 a) disks at ip — 0.79 
and b) diamonds at (p = 0.81. The trend of the 7* — 7** pairs is clearly linear; see Tab. Illll for parameters. The inset shows 7* 
vs Peo, with a guide to the eye (cubic polynomial in logPeo). 
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Table III. Parameters of the schematic model in the fits to the Brownian-dynamics simulation; the units are [t£] = fcsT '/d 2 a 
and [r] =D /d 2 a . 



overshoot) can be described. This was done by implementing a time-dependent prefactor (vertex) function in the 
generalized shear modulus git, 7), which is in accordance with the outcome of microscopic MCT as long as bare 
Peclet numbers are asymptotically small and the correlator decay is governed by strain. Th is concerns the pur e 
strain dependence of the prefactor function and its shape, which is motivated by the analysis bv lZausch et al. (2008). 
Effectively just one new parameter, namely a characteristic strain value 7*, is added to the model. It is the value of the 
strain where the maximal transient shear stress occurs. Thus, 7* can easily be read-off from experimental/ simulated 
data. It can also be calculated within the microscopic MCT, which however leads to values too large compared to the 
measured data. A second parameter of the schematic model, 7**, turned out to be linearly coupled to 7*. It describes 
the accumulated strain after which the stress relaxes onto the long-time limiting value, that makes up the flow curve. 

Brownian dynamics simulation data from a binary mixture of hard disks in two dimensions, and rheology measure- 
ments using thermosensitive core-shell microgels could quite well be analysed with the developed schematic model. 
Stress overshoots in the stress vs strain curves arose in the simulations and experiments even after eliminating waiting 
time dependences by equilibrating long enough. In the simulations, we could also verify that the flow remained ho- 
mogeneous throughout the start-up flow and in the steady state. The absence of these two (established) mechanisms 
being closely connected to non-monotonous stress-strain curves highlights that there is another intrinsic mechanism of 
the stress overshoot, which MCT identifies as negative stress auto-correlations. Microscopic MCT calculations reveal 
that structu ral distortion s with wave vectors corresponding to the average particle separation dominate during the 
stress peak ( Zausch et al.. 2008). Their anisotropy has not bee n clarified yet , but w as apparently contained to some 
degree in the isotropically sheared hard sphere model solved by IZausch et al\ ( 2008h . 

(Schematic) modeling of stress overshoots worked quite well with a single parameter set for a whole shear rate set, 
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albeit requiring us to adjust 7* as a weak function of shear rate. Comparison with microscopic MCT calculations in 
two dimensions suggest that the weak drift arises from a merging of the slow structural relaxation with the fast local 
motion. The model also misses a slow relaxation process which melts the glassy state. This process, often called the 
'hopping process' is missed by MCT in general, which predicts an ideal glassy state. As the stress overshoot arises 
in both the fluid and the glassy state, this deficiency of MCT, while noticeable here as well, does not rule out its 
applicability in the present context. 

As a future task, the close connection of stress overshoots and superdiffusive single particle motion, cf. IZausch et all 
(|2008h . can be studied within our schematic framework. We also plan to apply our stress overshoot description to 
measurement and simulation data of metallic glasses to check the universality of our generalized shear modulus form. 
Additionally, it will be highly interesting to investigate the features of the stress-strain curve together with the features 
of 7* and 7** starting from our microscopic MCT in three spatial dimensions. 
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